Monte Carlo simulation of pressure-induced phase transitions in spin-crossover 

materials 
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Pressure-induced phase transitions of spin-crossover materials were simulated by a Monte Carlo 
simulation in the constant pressure ensemble for the first time. Here, as the origin of the cooperative 
interaction, we adopt elastic interaction among the distortions of the lattice due to the difference of 
the molecular sizes in different spin states, i.e., the high spin (HS) state and the low spin (LS) state. 
We studied how the temperature dependence of the ordering process changes with the pressure, 
and we obtained a standard sequence of temperature dependences that has been found in changing 
other parameters such as strength of the ligand field (S. Miyashita et al.. Prog. Theor. Phys. 114, 
719 (2005)). Various effects of pressure on the spin-crossover ordering process are examined from a 
unified point of view. 

PACS numbers; 75.30.Wx, 74.62.Fj, 64.60.-1, 75.60.-d 
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Several spin-crossover (SC) compounds have been ex- 
tensively investigated 0,11, 11 S, H S 0, B S Ell , and 
various theoretical analyses of the SC transitions have 
been reported [H [H, [ll, [H, M, HI, H H E^. 

In the SC compounds, a metal ion can be in either 
a low-spin (LS) or high-spin (HS) state, depending on 
the strength of the ligand field. Control of the spin 
state of SC compounds has been realized by applying 
external stimuli such as temperature, light-irradiation 




i, 0, i S Jl|j_|i-, magnetic field 
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|22| . 1231 \2^ . and pressure (25. 26. 
has been pointed out that cooperative interactions play 
an important role for the SC transitions. With such inter- 
actions, various types of SC transition are realized, e.g., 
a smooth crossover or a discontinuous first-order phase 
transition. For modeling of the interaction mechanism, 
a model with an Ising-type short-range interaction, e.g., 
the Wajnflasz-Pick (WP) model, has been proposed, and 
various aspects of cooperative behavior have been sue- 



cessfuUy explained 



17, \li 



As the ligand field 



is changed, SC transitions show a sequence of tempera- 
ture dependences of HS fraction fasiT); (I) a smooth 
transition, (II) hysteresis, (HI) hysteresis with a low- 
temperature metastable HS phase, and (IV) a HS phase 
stable at all temperatures. We found that this sequence 
also appears with changing degeneracy or strength of the 
Thus we call this the generic se- 
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interaction 

quence. However, the origin of the interaction was not 
clear. Recently, it has been pointed out that the elastic 
interaction between distortions of the lattice due to the 
molecular size difference between the HS and LS states 
can induce phase transitions of the spin state [191] . In or- 
der to study the elastic interaction, besides the spin state, 
the positions of the molecules must be treated as dynam- 



ical variables to be equilibrated. This degree of freedom 
causes a change of the system volume. Therefore, the 
pressure P becomes an important parameter for describ- 
ing the system. Thus, we are now at a stage where we can 
study the pressure effect by direct numerical study. In a 
previous study [19], we demonstrated that the elastic in- 
teraction can cause a spin transition in a 2D system with 
an open boundary condition. In this study, we adopt a 
similar elastic model in 3D. Although intramolecular po- 
tentials were taken into account besides intermolecular 
potentials in Ref. ^1^, in the present study we take into 
account two levels (LS and HS) with two different molec- 
ular sizes as the molecular state for simplicity. Here we 
adopt a method for the constant pressure ensemble in 
a periodic boundary condition. Moreover, we adopt a 
Monte Carlo (MC) method where we can easily control 
the degeneracy of the spin state. By this method, we 
succeeded for the first time in demonstrating pressure- 
induced phase transitions in a spin-crossover system. 

The pressure effect has been one of the most important 
characteristics observed in SC compounds, and various 
interesting properties have been observed, e.g., a shift of 
the transition temperature, increase or decrease of hys- 
teresis, stabilization of the LS state for the whole tem- 
perature region, and so on [25, 26, 27, 28, 29, 30|]. To un- 
derstand the pressure effects, a theoretical analysis based 
on each experimental phenomenon and the comprehen- 
sive theory of the free energy of the mean-field model has 
been reported H, 13, 1^. In the present work, we study 
these effects by a direct numerical method, adopting a 
microscopic Hamiltonian and demonstrating the funda- 
mental aspects in a unified microscopic picture. 

We performed MC simulation on the simple cubic lat- 
tice. In the spin-crossover materials, the molecule at a 
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LS molecule 




Spring constant k; 



FIG. 1; (Color online) Schematic illustration of the present 
model. HS/LS molecule consists of Fe atom (red/blue circle) 
and six ligands (gray circles). 



lattice point may be in the HS or LS state. We express 
the spin state at the z-th site by n.i, which equals for 
the LS state and 1 for the HS state. As important ingre- 
dients of the spin-crossover material, we set the energy 
difference of the state D{> 0), and the degeneracies of 
the states: gns and gLS for the HS and LS states, respec- 
tively. These properties are represented by the on-site 
Hamiltonian 



rii, 



(1) 



In order to introduce interaction between spins, an Ising- 
like interaction —JJ2<ij> f^if^j was adopted in the WP 
model. Instead, we have attributed the interaction to 
the elastic interaction between lattice distortions caused 
by the difference of the molecular size between the HS or 
LS states. Therefore, we introduce the elastic interaction 
between the molecules: 



nnn i' -^-i-O, 
2 



(2) 
(3) 

(4) 



where is the distance between the molecule on the i- 
th and j-th sites, and ki is the corresponding spring con- 
stant (Fig. [1]). Hnn expresses elastic interaction between 
nearest-neighbor pairs The interaction is a func- 

tion of the molecular radius Ri{= (1 — ni)Rhs + riiR^s) 
of the i-th site, where i?HS and RhS are the molecular ra- 
dius of HS and LS states, respectively. We set the ratio 
of the radii to be i?Hs/-RLS = 1.1- i?nnn expresses elastic 
interaction for next-nearest neighbor pairs In 
this study, we set the ratio of the spring constants, fci/fc2, 
to be 10 

For the simulation, we adopt the NPT-MG method 
1321 for the isothermal-isobaric ensemble with the number 



of molecules N, the pressure of the system P, and the 
temperature T. The thermodynamic potential for the 
isothermal-isobaric ensemble is the enthalpy. 



H = E + PV, 



(5) 



where E is the energy given by Eq. (1) and V is the vol- 
ume of the system. The states of the system are specified 
by AN + 1 variables (ni, • • • , njv, ri, • • • , rjv, F). In the 
NPT-MC method, we have the following balance condi- 
tion for the transition probabilities {wi^k}- 



Wk- 



exp(-/3AVF), 



(6) 



where 



AVF = {Ek - E,) + P{Vk - V-) - iVriog(-^). (7) 

The scheme of the simulation is as follows: (i) Choose 
a molecule randomly, (ii) Choose a candidate spin state 
m ^ or 1 by the probability 5ls/(5ls + 5hs) or 
5hs/(3ls + 5hs), respectively, (iii) Choose a candidate 
position of the molecule (r^). Here, we use the scaled 
coordination length s = r^/i {L ~ V^^^). We choose the 
candidate position as 



Sz + S^,, (8) 



where S^xi S,y and are random numbers between —1 and 
1, and 5 0.005. (iv) Update the state by the Metropo- 
lis method, (v) Repeat the above update N times (vi) 
Choose a candidate for a new size of the system with a 
random number 



(9) 



(vii) Update the size L. Here, A is taken equal to 
O.OSTVi/^i^Ls. In this study, we performed 10000 MCSs 
for transient steps and 10000 MCSs to measure the physi- 
cal quantities. The system size is iV = 16 x 16 x 16, which 
is enough to study thermal properties. 

First, we study how the types of temperature depen- 
dence of the HS fraction /hsC^) change with the spring 
constant fci(= 10fc2). In Fig.[2]we depict f-as{T) for var- 
ious values of fci with g{= 5hs/5ls) = 20, D = 1,P = 
0.01, and fcl = 10,20,30,40, and 50. When ki is small, 
e.g., fci = 10, the transition is gradual and hysteresis 
is not observed. As fci becomes large, the transition be- 
comes sharp. For fci = 40, hysteresis is observed between 
cooling and warming processes. Here, the transition tem- 
perature is 0.28 for the cooling process and 0.36 for the 
warming process. Snapshots of the spin configuration 
(a two-dimensional section of the lattice) are shown in 
Fig. [3] When fci becomes larger, i.e., fci = 50, the tran- 
sition does not take place in the cooling process and the 
HS phase is maintained down to T = 0. This observa- 
tion indicates the existence of a HS metastable phase at 
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FIG. 2: (Color online) Temperature dependence of the HS 
fraction /Hs(r) with g = 20, L> = 1,P = 0.01, and ki = 10 
(red squares), 20 (green circles), 30 (light blue diamonds), 40 
(blue triangles), and 50 (pink inverted triangles). Solid lines 
are guides for the eye. 



FIG. 4: Temperature dependence of the HS fraction with 
g = 20, D = 1, P = 0.01, and fci = 40. Dotted line is a guide 
for the eye. 
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T = 0.16 
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FIG. 3: (Color online) Two-dimensional snapshots at various 
temperatures with g = 20, D — 1,P — 0.01, and fci = 40. 
Gray (red online) and black (blue online) circles are HS and 
LS molecules, respectively. 



low temperatures [17|, . In the warming process from 
T — 0, the LS phase is transformed to the HS phase at 
T = 0.4. 

We find that the change of fci causes a sequence of 
fusiT) that agrees with the generic sequence proposed 
in our previous papers. Thus, we expect that there is a 
case where a low-temperature HS metastable phase and 
thermal hysteresis are observed, which was found charac- 
teristic of this type of ordering processes and also was ex- 
perimentally confirmed IT], llSl . l33| . In order to check the 
existence of the low-temperature HS metastable phase, 
we studied the warming-up process from HS from T = 
for the system with fci = 40. The temperature depen- 
dence is depicted in Fig. [H where the HS metastable 
phase exists and relaxes to the LS phase at T = 0.04. 
The LS phase changes to the HS phase at T = 0.36 as we 
saw in Fig. [51 Therefore, we find that the HS metastable 
phase and hysteresis are both observed. Now, we con- 
firm that the present model realizes the generic sequence 
of /hs(7") of the SC transitions. 

Next, we study the pressure effect on the SC transi- 
tions. We study how /hs(2^) for fci = 50 changes when 
the pressure increases. In Fig. [51 we depict fusiT) dur- 




FIG. 5: (Color online) Temperature dependence of HS frac- 
tion with g — 20, D = 1, fci = 50, and P — 0.01 (pink inverted 
triangles), 0.1 (red squares), 0.2 (green circles), 0.5 (blue tri- 
angles). Inset is the pressure dependence of T1/2. Solid lines 
are guides for the eye. 



ing the warming process from the HS state. In the case 
of low pressure, P — 0.01, the transition is not observed, 
as shown in Fig. [2l When the pressure becomes P = 0.1, 
the initial HS state relaxes to the LS phase at T = 0.04, 
which indicates a low-temperature metastable HS state. 
Then, we find hysteresis with a jump at T = 0.32 in the 
cooling process and at T = 0.44 in the warming process. 
For P = 0.2, tlie initial HS state immediately relaxes 
to the LS phase, which indicates no low-temperature 
metastable HS state. Here, the hysteresis disappears and 
fns{T) is shifted to the high-temperature side. When 
the pressure increases further high, P = 0.5, the transi- 
tion becomes gradual, and fnsiT) is shifted to tlie high- 
temperature side. We define T1/2 at which /hs — 1/2. 
In the hysteresis region we define it as (1/2 — /hs) in the 
warming process is equal the (/hs ~ 1/2) in the cooling 
process. The pressure dependence of T1/2 is depicted in 
the inset of Fig. [5l which indicates that ri/2 increases 
linearly with the pressure. 

We also study the pressure dependences of /hs (T) at 
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Pressure 

FIG. 6; (Color online) Pressure dependence of the HS frac- 
tion with g = 20, D = 1, ki = 50, and T — 0.3 (red squares), 
0.5 (green circles), and 0.7 (blue triangles). The numbered 
arrows indicate the pressure sequence. Solid lines are guides 
for the eye. 

various temperatures. In Fig. [6l fusiT) a.t T ~ 0.3,0.5, 
and 0.7 for fci — 50 are depicted as functions of the pres- 
sure. Here, the HS phase is set as the initial phase. At 
T = 0.3, a small pressure P = 0.1 induces the transi- 
tion from the HS to the LS phases (red (i) arrow). The 
pressure-induced LS phase does not return to the HS 
phase in the process of reducing pressure (red (ii) ar- 
row). This irreversible pressure effect indicates that the 
pressure stabilizes the LS phase and destabilizes the HS 
phase. At T = 0.5, the transition from the HS to the 
LS phase is observed at P = 0.22 during the pressure- 
increasing process (green (i) arrow). During the pressure- 
reducing process, the transition from the LS to the HS 
phases is observed at P = 0.16 (green (ii) arrow). That 
is, in this case we observe pressure-induced hysteresis. At 
T = 0.7, the transition between the HS and LS phases is 
smooth, and the hysteresis disappears. The present ob- 
servations indicate that the pressure plays a similar role 
to that of the temperature for the SC transition. 

The NPT-MC method for a three-dimensional system 
was established, and the effect of the pressure on the SC 
transition was studied. To our knowledge, this is the first 
attempt to study the pressure effect by direct numerical 
simulation, considering the local lattice distortions cased 
by the molecular size difference between HS and LS states 
in SC complex. In particular, we succeeded in observing 
a sequence of /hs (T) as a function of the pressure which 
agrees with that proposed as a general sequence the SC 
transitions. 

In the present study, we kept the parameters D and 
g constant. If we take into account this pressure depen- 
dence, we can have a great variety of pressure depen- 
dences, which correspond to complicated dependences 
observed in experiments. For the next stage, we will 
study various pressure effects from the viewpoint of the 
present model and attempt to obtain a systematic un- 
derstanding of the variety of pressure effects on the SC 



transitions. 
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